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Abstract 

We use density functional theory (DFT) with the generalized gradient approximation (GGA) 
and our first-principles extrapolation method for accurate chemisorption energies [Mason et al, 
Phys. Rev. B 69, 161401R (2004)] to calculate the chemisorption energy for CO on a variety 
of transition metal surfaces for various adsorbate densities and patterns. We identify adsorbate 
through-space repulsion, bonding competition, and substrate-mediated electron derealization as 
key factors determining preferred chemisorption patterns for different metal surfaces and adsorbate 
coverages. We discuss how the balance of these interactions, along with the inherent adsorption 
site preference on each metal surface, can explain observed CO adsorbate patterns at different 
coverages. 
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I. INTRODUCTION 



The prototypical CO/transition metal adsorption systems, though extensively studied, 
continue to be of significant contemporary research interest. These systems owe 

their continued examination in part to the significant variations in CO bonding properties 
from one transition metal to another, such as the evolution of chemisorption geometry 
with adsorbate coverage 0. Attempts to describe adsorption trends as a function of are 
complicated by the many phenomena involved, such as through-space interactions between 
adsorbates, substrate-mediated interactions, surface relaxations and reconstructions induced 
by adsorption, adsorption at different high-symmetry sites, patterned overlayer formation, 
and transitions between different adsorbate arrangements on the surface under different 
conditions. 



There are many studies that address adsorbate-adsorbate interactions 
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14 ll^ , and it is beyond the present scope to review these and other studies in full. 
One major conclusion is that properties of sub-monolayer CO overlayers such as chemisorp- 
tion energy and surface site occupation are found to be significantly coverage dependent. 
Interactions between CO molecules are strongly repulsive at short range; beyond nearest 
neighbor distance, inter-adsorbate interactions are non-monotonic and vary as a function 
of adsorbate coverage. A better understanding of adsorbate interactions would be helpful 
for understanding the wide va riety of adsorption systems of interest to the surface science 
community. Q fl 3, Q, Q Q 

To understand how chemisorption evolves with 0, it is necessary to identify the factors 
influencing the adsorbate-adsorbate interactions, and how these factors vary from metal to 
metal. For this purpose, it is advantageous to have a data set spanning a variety of metals, 
and with a broad range (in terms of coverage, site occupancy, and overlayer pattern) of 
chemisorption geometries. With such a dataset, chemisorption energies between different CO 
overlayer structures can be compared for different metals, elucidating trends and revealing 
important interactions. 

The goal of comparing chemisorption at different coverages, and at different adsorbate 
overlayer patterns for a given coverage, is well suited for theoretical study. We can model 
adsorption at coverages and in patterns not experimentally observed. However, it is nec- 
essary to have a means for accurately modeling chemisorption on a variety of metals, and 
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for all adsorption sites, coverages, and overlayers of interest. Owing to advances in com- 
puter power and algorithms, many aspects of surface science can be quantitatively explained 
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through theoretical studies at the DFT-GGA level. 

Hindering the use of DFT for quantitative investigation of CO/metal adsorption systems 
is the tendency of DFT calculations to overestimate chemisorption energies and to favor 
multiply coordinated adsorption sites over the top site, in disagreement with experimental 
data. This "CO/metal puzzle" was first comprehensively addressed by Feibelman et 0,1^2^ 
Subsequently, several studies have traced the CO/metal puzzle to its origin. [24. 125L 126 1 

Following these investigations, we recently derived a first-principles extrapolation pro- 
cedure 0] which, when applied to eight different single-crystal substrates, reduced over- 
binding and achieved site preferences consistent with experimentally observed adsorption in 
all cases. Other works discussing how to correct for the DFT site-preference problem include 
Olsen et al. 28] and Kresse et al. 2(| In the work of Olsen et al, relativistic effects are found 
to be important to correct site preferences for CO on Pt (111). However, this fails to explain 
why DFT finds incorrect site preference on non-relativistic Cu, as recently also pointed out 
by Gajdos et al. [2^ Kresse and co-workers employ a GGA+U type exchange-correlation 
energy functional. However, the correct value of U is not known a priori. In contrast, our 
method uses ab initio calculations to unambiguously calculate the correction magnitude. 
The present work is the first time any of these three accurate methods has been applied to 
multiple metals, surfaces, and coverages in a single study. 



II. METHODOLOGY AND COMPUTATIONAL DETAILS 

In the present study, we investigate chemisorption of CO on transition metals as a function 
of metal, surface termination, adsorption site, and CO overlayer pattern. The focus of 
the present work is transition metal chemisorption, but some additional results for CO on 
Al (111) and (100) surfaces are also presented for the purpose of contrasting transition 
metals with an sp metal. 

DFT calculations are carried out using the PBE GGA exchange-correlation functional 
and norm- conserving optimized pseudopotentials J^l] with the designed nonlocal method 
for metals. 3, [3] All pseudopotentials were designed using the OPIUM pseudopotential 
package. |24| Metal surfaces are modeled as c(4 x 2) or p{2 x 2) slabs of five layers separated 
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by vacuum, with atomic relaxation allowed in the two topmost layers. Calculations are 
done, and E^em tested to be converged within less than 0.025 eV, using a 4 x 4 x 1 grid of 
irreducible Monkhorst-Pack fc-points. [3^ 

All values for E L ^ m reported here have been corrected using our first-principles extrapo- 
lation procedure. 27[ We have tested that the corrections to E chem do not appreciably evolve 
with 9, affecting corrected values for E c ^ era by less than 0.01 eV/molecule. CO chemisorp- 
tion is modeled at coverages of 0=0.25, O=0.5, 0=0.75, and 0=1 ML (monolayer). 

In these calculations, the CO bond is held perpendicular to the surface. Eliminating 
the CO tilt angle as a degree of freedom facilitates comparison of chemisorbed and gas- 
phase adsorbate-adsorbate interactions, as well as comparison of CO-CO interactions over 
different metals. None of the higher-coverage chemisorption structures on the (100) surfaces 
were found to have lateral forces on CO under this constraint. In some of the higher-coverage 
structures on the (111) surfaces, lateral forces of psO.01-0.05 eV/A were found on the C and 
O atoms, as indicated in Table HI No dipole correction scheme to decouple the periodic slabs 
was included, as such corrections were found to affect the total energy of the chemisorption 
systems by less than 0.01 eV/cell. 

Values for -E^em are calculated as 

E c hem = (— E 

surface— CO 

+ NE co + E suriace )/N (1) 

where N is the number of CO molecules per unit cell. With this definition, chemisorption 
energies are reported in eV/CO molecule, and positive values of E chern represent favorable 
adsorption on the surface. 

Note: The overlayer naming convention denotes both coverage and sites: one letter for 
0=0.25, two for O=0.5, three for 0=0.75, and four letters for 0=1. 

On the (111) surfaces of Pt, Rh, and Pd, we modeled chemisorption at 0=0.25 at top 
(t), bridge (b), hep (h) and fee (f) adsorption sites, as shown in Figure [1(a)] We consider 
our results for 0=0.25 to be the low-coverage limit values of -E c hem- As a test of the validity 
of this assumption, we also calculated -E c hem o n each metal at each site in a p(2 x 3) cell, 
corresponding to a coverage of 0=0.17, on both the (111) and (100) surfaces. The change 
in i?chem from 0=0.17 to 0=0.25 was less than 0.04 eV for all sites and metals, with the 
same site preference. The most important characteristic of the t, b, f, and h arrangements 
is the absence of nearest neighbors. 
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At 0=0.5, we modeled chemisorption at top (tt) and hep (hh) adsorption sites, as shown 
in Figure |T(b)[ We did not model fee adsorption at 0=0.5 due to the very similar energies for 
f and h for all surfaces as shown in Table The tt and hh structures allow us to investigate 
how chemisorption changes when adsorbates are separated by nearest neighbor distances. 

We model three different bridge site configurations, bbl, bb2, and bb3, as shown in 
Figures |l (b)[ |1 (c)| and |l(d)| In bbl and bb2, nearest neighbor CO molecules are present. For 
the bb2 configuration there is an additional characteristic: each surface metal atom involved 
in chemisorption participates in two carbon-metal bonds, leading to bonding competition, |4| 
The bb3 configuration does not have nearest neighbor CO molecules, and is one of the 
observed patterns at 0=0.5 on the Pd (111) surface. 3(| Comparing these three bridge 
patterns allows us to deduce the energetic costs of nearest-neighbor adsorbate interactions 
and bonding competition for each metal. 

We also considered (see Figure |2(a)|) mixed occupation of top and bridge (tb) in a c(4 x 2) 

nnn 

cell, the observed pattern at 0=0.5 on Pt (111) [23, |3ja, |3J|, and mixed occupation (see 
Figure |2(bj| of hep and fee (hf) in ap(2x2) cell, one of the observed patterns on Pd (111). |36l | 
For the hf over layer, bonding competition between CO molecules is present. 



At 0=0.75, we modeled an overlayer with CO in top, hep hollow, and fee hollow 
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Figure p(b)[ This coverage and pattern is observed on both Rh and Pd (111) 
and is another configuration where bonding competition is present. 

Chemisorption at 0=1 is modeled in a c(4 x 2) cell for both top site (tttt) and hep (hhhh) 
adsorption. Bonding competition cannot exist between top site adsorbed molecules, but is 
present in the hhhh configuration. Consistent with previous work, we find chemisorption at 
0=1 to be exothermic. The lack of experimental observation of CO at 0=1 is due to kinetic 



effects. 



42] 



On the (100) surfaces of Pt, Rh, and Pd, we modeled chemisorption at 0=0.25 at top 
(t), bridge (b), and hollow (h) adsorption sites (see Figure |3(a)| . As with the (111) surfaces, 
we consider adsorbates at this coverage to be isolated. 

At 0=0.5 on (100) surfaces, we modeled chemisorption at two different top site con- 
figurations, ttl and tt2, two different hollow site configurations, hhl and hh2, and three 
different bridge site configurations, bbl, bb2, and bb3. These overlayer patterns are shown 
in Figures |3(b)||3(c]j andgd)) The top sites vary in that ttl has CO molecules at nearest 
neighbor distances, while tt2 does not. Likewise, the hhl overlayer has adsorbates at nearest 
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(a) 



(b) 




(c) (d) 

FIG. 1: (a) Schematic of (111) surfaces with CO at 0=0.25. Overlayer t is indicated by squares, b 
by circles, and h by "X". The cell is indicated by the shaded region, (b) Schematic of (111) surface 
with CO at 0=0.5 in overlayer tt (squares), in the bbl overlayer pattern (circles) and hh pattern 
("X"). (c) Schematic of the (111) bb2 structure, (d) Schematic of (111) bb3 structure. 




(a) (b) 

FIG. 2: (a) Schematic of (111) tb structure at 0=0.5. Occupied top sites are shown with squares, 
and occupied bridge sites are shown with circles, (b) Schematic of (111) hf and thf structures at 
O=0.5 and 0=0.75. In the former, only hep and fee sites ("X") are occupied. In the latter, the 
top sites (squares) are also occupied. 
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FIG. 3: (a) Schematic of (100) surface with CO at 0=0.25. Top site pattern t is indicated by 
squares, overlayer b by circles, and h by "X". (b) Schematic of (100) ttl, bbl, and hhl overlayers 
(6=0.5). (c) Schematic of the (100) tt2, bb2, and hh2 overlayers. (d) Schematic of (100) bb3 
structure with CO at 0=0.5. 



neighbor separations while the hh2 overlayer does not. In the bridge site structures, both 
bbl and bb2 have CO molecules at nearest neighbor separations, but differ in that bb2 has 
bonding competition while bbl does not. The bb3 pattern is experimentally observed on 



both Pt and Rh (100). |43l |44j We also modeled tttt and hhhh adsorption at 0=1. All 



modeling of CO on the (100) surfaces was done in a c(4 x 2) cell. 



III. RESULTS 



We begin with low coverage (0=0.25) results. Detailed comparison of these calculated 
energies and experimental results have been reported previously. [2^ In the low-coverage 
limit, our calculations reproduce the site preferences observed experimentally for Rh, Pd 
and Pt (111) surfaces (top for Rh and Pt, fee for Pd). 15] Experimental data show that 
the lowest-coverage ordered CO structure on all three of these (111) surfaces is (^/3 x 
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V / 3)-R30° at 0=0.33. The formation of this overlayer is driven by nearest neighbor CO-CO 
repulsive interactions, since this overlayer gives the highest- coverage single-site occupation 
without nearest neighbor interactions. The site occupied in the (\/3 x y/3)R30° structure is 
determined by the site preference of isolated CO adsorption. 

The (100) surfaces of Pt, Rh, and Pd share common low-coverage chemisorption behavior. 
The lowest coverage at which ordered CO structures are observed on the (100) surfaces is 
O=0.5. For (100) surfaces, we find that difference in E c h em between the 0=0.25 CO overlayer 
and the highest E c ^ em O=0.5 CO overlayer is negligible, in agreement with experiments that 
find a c(4 x 2) overlayer at O=0.5 to be the lowest coverage ordered structure. Q] On 
Pt (100), domains of top-site CO and of bridge-site CO are seen. On Rh (100), CO 
occupies top sites at O=0.5, and on Pd (100) CO binds at bridge sites. [3] For all three 
metals, our results for site preference at 0=0.25 are consistent with these experimental 
observations. 

While Pd, Rh and Pt (111) surfaces form the same ordered structure at low coverage, 
they evolve differently with coverage. Comparison between DFT and experiment at higher 
O is imprecise due to our neglect of temperature, pressure, and entropy effects. However, our 
DFT results for O=0.5 and K are consistent with experimental results. Our calculations 
show that chemisorption of CO on Pt (111) at O=0.5 in the c(4 x 2) tb mixed occupation 
structure is favored by « 0.09 eV/molecule over the next highest -E c hem overlayer, i n ag ree- 
ment with the experimental observations of the c(4 x 2) tb structure at O=0.5. p\. ws. I39I] 
For Pd (111) we find that the bb3 c(4 x 2) structure is favored by « 0.04 eV over the c(4 x 2) 
hf structure. All other overlayer patterns show -E c hem values more than 0.1 eV lower than 
the hf and bb3 structures. This result is in agreement with coexistence of c(4 x 2) bb3 and 

n 

hf overlayers at O=0.5 observed on Pd (111) in low-temperature STM experiments. |3fi| 
The case of Rh (111) is especially interesting, with the experimental LEED pattern 



showing disorder at O=0.5. 



4l| Our results show that the bb3, tb, and hf overlayers on 



Rh (111) at O=0.5 have the same -E c hem to within 0.03 eV. The experimental result may be 
due to full disorder, but our computations lead us to suggest that domains of several nearly 
isoenergetic overlayer patterns coexist, explaining the complicated LEED pattern. 

For Pt (100) and Rh (100) surfaces, we find a small energy difference between the top and 



bridge site, consistent with the occupation of both of these sites for O=0.5 and greater. 



43, 



1441 ] For Pd (100), the large site preference energy is consistent with experimental observation 
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TABLE I: Results for CO E chem at different O on the (111) surfaces of Pt, Rh, and Pd, in eV. 
Selected Al values are also reported. For coverages above 0=0.25, the difference E- mt between 
the low-coverage value and high-coverage value -Echem for the corresponding sites is reported in 
parenthesis. Values for experimentally observed structures are in italics. * indicates systems for 
which the lateral forces on CO were not minimized (CO is constrained to be perpendicular in all 
cases) . 





Pt (111) 


Rh (111) 


Pd (111) 


Al (111) 




-^chem -^int 


-^chem -Eint 


-^chem -^int 


-^chem -^int 


0=0.25 










t 


1.56 


1.67 


1.25 


0.22 


1) 


1.43 


1.58 


1.49 




h 


1.40 


1.64 


1.60 




f 


1.43 


1.64 


1.63 




0=0.5 










tt 


1.27 0.29 


1.35 0.32 


0.93 0.32 


0.17 0.05 


bbl* 


1.12 0.31 


1.25 0.33 


1.15 0.34 




bb2 


1.07 0.36 


1.40 0.18 


1.09 0.40 




bb3* 


1.36 0.07 


1.55 0.03 


1.44 0.05 




tb* 


145 0.04 


1.56 0.06 


1.31 0.06 




hh* 


1.03 0.37 


1.33 0.31 


1.19 0.41 




hf 


1.21 0.20 


1.53 0.11 


1.40 0.21 




0=0.75 










thf* 


1.19 0.27 


1.39 0.26 


1.15 0.34 




0=1 










tttt 


0.74 0.82 


0.53 1.14 


0.39 0.86 




hhhh* 


0.31 1.09 


0.76 0.88 


0.46 1.14 





that only bridge sites are occupied for all values of 0. |4 
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TABLE II: Same as Tabled for CO chemisorption on (100) surfaces. 





Pt (100) 


Rh (100) 


Pd (100) 


Al (100) 




-^chem -^int 


-^chem -^int 


-^chem -^int 


-^chem -^int 


G=0.25 










t 


1.80 


1.72 


1.35 


0.12 


b 


1.83 


1.77 


1.64 




h 


1.25 


1.60 


1.50 




6=0.5 










ttl 


1.49 0.31 


1.39 0.33 


1.01 0.34 


0.07 0.05 


tt2 


1.78 0.02 


1.72 0.00 


1.34 0.01 




bbl 


1.60 0.23 


1.47 0.30 


1.36 0.28 




bb2 


1.50 0.33 


1.58 0.19 


1.28 0.36 




bb3 


1. 77 0.06 


1.73 0.04 


1.62 0.02 




hhl 


0.83 0.42 


1.23 0.37 


1.10 0.40 




hh2 


1.08 0.17 


1.49 0.11 


1.35 0.15 




6=1 










tttt 


1.16 0.64 


1.08 0.64 


0.70 0.65 




hhhh 


0.11 1.14 


0.73 0.87 


0.38 1.12 





IV. DISCUSSION 



In this section, we discuss trends in the data presented in Tables |U and |TT| We identify 
the presence of nearest neighbors and/or metal-atom sharing as the chief determinants of 
interactions between adsorbed CO molecules. 

The effects of the substrate on adsorbate-adsorbate interactions will be discussed in the 
framework of the Hammer-Morikawa-N0rskov model 4|| for the <i-band contribution to 
top-site CO chemisorption: 



+2 



(1-/)- 



^271-* — 

V? 



+ fs n v n 



+ (1 + f)S a V a 



(2) 
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where / is the idealized filling of the metal d bands, V% is the metal-carbon orbital overlap 
matrix element at the top site, €2n* and are the energies of the renormalized 2n* and 5a 
orbitals relative to the Fermi level, and is the center of the metal <i-bands relative to the 
Fermi level. V and S are coupling and overlap matrix elements respectively, and are labeled 
by orbital. 

We also analyze gas-phase CO-CO interactions, study induced charge densities in 
chemisorbed systems, and analyze the metal electronic structure. 

A. Nearest-neighbor interactions and chemisorption-induced perturbations 

We now analyze the E c ^ em results presented in Section III to understand nearest-neighbor 
adsorbate-adsorbate interactions, for systems that have no metal-atom sharing. 

We can estimate the energetic cost of placing adsorbates at nearest neighbor distances 
from our top-site chemisorption results. Our dataset includes three different exclusively 
top site chemisorbed structures at 0=0.5: The tt structure on (111) (see Figure |l(b)| ), 
and the ttl and tt2 structures on the (100) surface (see Figures ggj] and gcj . While the 
chemisorbed CO molecules in the tt and ttl structures have nearest-neighbor CO adsorbates, 
the adsorbates in the tt2 structure do not. From inspection of the data in Tables HI and HTl it 
is clear that E int (the difference between the low-coverage limit and higher coverage -E c hem) 
for the tt2 structure on the (100) surface is tiny (<0.02 eV) for all metals. E int is much larger 
but still nearly independent of metal identity for the tt and ttl patterns (=0.32 eV±0.02 eV). 
E int for the Al surfaces is almost an order of magnitude smaller than for the transition metal 
surfaces. 

We can compare these E int results with the gas-phase CO-CO interaction energy Ef^, 
calculated by: 

ES:(r) = E^_ co -2E co (3) 

The interaction energy for CO at different separations is plotted in Figure 0] At the 
nearest-neighbor separations for Al (2.89 A), Pt (2.76 A), Pd (2.73 A), and Rh (2.66 A), the 
gas-phase CO-CO interaction energies are 0.12, 0.19, 0.21, and 0.27 eV, respectively. The 
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FIG. 4: Calculated interaction energy Ef^(r) for 2 CO molecules at different separations r. 

corresponding (111) tt E- mt values from Table|I]are 0.05, 0.30, 0.32, and 0.32 eV, respectively. 
For Al, E int < Ef^, suggesting that the s and p metal electrons do not increase CO-CO 
repulsion and perhaps lessen it somewhat. The transition metals all show Ei nt > Ef^, 
meaning that through-space repulsion is augmented by the surface. From the contrast 
with Al, we deduce that this increase in repulsion must be due to d electrons changing the 
adsorbed CO electronic structure. 

The Hammer- Morikawa-N0rskov model (Equation^ shows that the back-donation from 
the metal <i-bands into the CO 2n* increases with /. Increased filling of the 2it* or- 
bitals strengthens the repulsion between chemisorbed CO molecules, while increased inter- 
adsorbate separation weakens it. Pt and Pd have a higher / but a greater CO-CO separation 
than on Rh. Our data indicate that these effects balance, giving nearly the same CO-CO 
interaction energy (top-site E int ) for these three transition metals. 

The bridge-site systems show the same E int trends. For systems with no first nearest 
neighbors (bb3) E{ nt is very small. For bridge overlayers at O=0.5 with nearest neighbors 
but no metal-atom sharing (bbl), a large metal-independent E int w 0.32 eV is observed, 
just as for top sites. (Actually, the bbl pattern on Pt (100) is somewhat more stable; this 
surface is known to reconstruct, so chemisorption may be stabilized by relieving surface 
48,|3] The changes in E chem we do observe on the unreconstructed Pt(100) 



stress. 43 
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surface are in excellent agreement with other recent theoretical work. |3j) 

Nearest-neighbor hollow sites always involve metal atom sharing, so they are considered 
in the next section. 
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TABLE III: Values for the d-band center, e<£, for the (111) and (100) surfaces of Pt, Rh, and Pd 
when no CO is adsorbed (0=0), and when CO is adsorbed at bridge site at 0=0.25. The changes 
in Ed due to chemisorption (A) are also listed. All reported values are in eV. 





O=0 


=0.25 


A 


Pt (111) 


-2.93 - 


-3.59 


-0.66 


Pt (100) 


-2.79 ■ 


■3.50 


-0.71 


Pd (111) 


-2.29 ■ 


■3.04 


-0.75 


Pd (100) 


-2.17 ■ 


■2.94 


-0.77 


Rh (111) 


-2.42 - 


-2.85 


-0.43 


Rh (100) 


-2.29 - 


-2.73 


-0.44 



B. Metal-atom sharing: bonding competition and electron derealization 

Bonding competition arises when the surface metal atoms involved in chemisorption par- 
ticipate in more than one carbon-metal bond. To find the susceptibility of -E c hem to external 
perturbation (such as the sharing of metal atoms) within the Hammer-Morikawa-N0rskov 
model, we differentiate the expression for Ef hem (Equation |2J) with respect to external pa- 
rameter A, representing the number of chemisorption bonds per metal atom: 



dE. 



chem 



dX 



2(1 -m 



de d 
dX 



(4) 



[e 27r * — td) \t d — e 5cr j 

The opposite signs of the two terms in Equation 0] mean that any change in the position 
of €d has two competing effects on E chern . This competition vanishes as / — > 1, so E chem on 
Pd or Pt (/ =0.9) should be more sensitive to bonding competition than Rh (/ =0.8). 

The change of -E c hem with bonding competition A also depends on ded/dX, the shift in the 



rf-band energy as chemisorption bonds are 


"ormed. The e, 


i position is affected by perturbation 


to the metal surface such as strain 


50, 


51 


|, alloying [ 


3 


m 


] and creation of chemical bonds 



with adsorbates. 4j, |53( This effect is explored in Table 11111 showing that CO chemisorption 
causes a much smaller shift in e d for Rh than for either Pd or Pt. 

The preceding analysis of the Hammer- Morikawa-N0rskov model suggests that bonding 
competition should increase E- m % for all chemisorption involving metal sharing, but much 
less on Rh than on Pt or Pd, due to lower / and lower ded/dX. 
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<f> < j > 



(a) 




FIG. 5: Charge density changes induced by CO adsorption in pattern bbl site on Pt (100). Lightly 
shaded iso-surfaces show areas of charge gain, and darker iso-surfaces show areas of charge loss, 
(a) Top view, (b) Side view. 

In addition to the destabilizing effect of bonding competition, metal-atom sharing also 
enables electron derealization. As illustrated in Figures and El the CO 2tt* orbitals of 
neighboring molecules and the metal d orbitals of the atom they share combine to form 
extended states in the bb2 system (Figure EJ). This derealization is not present on the bbl 
overlayer (Figure EJ) • 

For the bridge-bonded systems with metal-atom sharing (bb2), the effect is dramatic. 
On Rh (111) and (100), derealization stabilizes chemisorption significantly, and bonding 
competition gives a only small reduction, leading to stronger chemisorption and E [nt values 
that are 0.1-0.15 eV lower for bb2 than for bbl. On Pt and Pd, much stronger bonding 
competition destabilizes chemisorption, raising E- int for bb2 about 0.15-0.2 eV above that of 
Rh, and 0.05-0.1 eV above the bbl values on Pt and Pd despite electron derealization. 

Hollow site chemisorption exhibits the same trends, with bonding competition mak- 
ing bonding less effective, and metal-mediated derealization stabilizing metal-sharing 



14 



CI $ I 




(b) 

FIG. 6: Same as Figure [51 for the bb2 overlayer. Note the significant charge sharing between 
neighboring CO molecules, mediated by the shared metal atom. 

chemisorption. However, because a CO at a hollow site distributes its bonding among 
three or four surface atoms, both the bonding competition and derealization effects from 
metal sharing are somewhat less than for bridge sites. In all cases (hh,hhl,hh2) E- mt for Rh 
is more favorable by ~0.05 eV than on Pt and Pd. Since E int at hh and hhl is greater than 
at tt and ttl, metal sharing appears to always be net repulsive for hollow sites. 



C. Preferred Overlayer Structures 

In this section we discuss how adsorbate-adsorbate interactions, along with the low- 
coverage CO site preference on each surface, influence -E^em for various 0=0.5 overlayers. 

On the (111) surfaces of Pt, Pd, and Rh, we modeled chemisorption at 0=0.5 for seven 
different ove rlayer p atterns. In the case of Pt (111), we find tb (experimentally observed 

on pt (in) 0HQ) 

has the highest -Echem, by at least 0.09 eV over the other overlayer 
patterns. At low coverage, t is preferred by 0.13 eV over b and f, and by 0.16 eV over h. 
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On the (111) surface, exclusive occupation of top sites (at >0.5) must involve nearest- 
neighbor adsorbates, whereas the mixed occupation tb avoids nearest neighbors, reducing 
through-space repulsion. Therefore, the tb structure can be viewed as a compromise between 
the site preference for t, and reduced adsorbate-adsorbate repulsion made possible by partial 
b occupation. 

n 

On Pd (111), we find -Echem for bb3 and hf (both of which are experimentally observed 36]) 
to be within 0.04 eV of each other, and at least 0.09 eV stronger than the other five patterns. 
At low coverage, we find f to be preferred by 0.14 eV over b. Here we can also interpret 
the hf and bb3 overlayer patterns to be compromises between inherent site preference and 
adsorbate-adsorbate repulsion minimization. In the bb3 pattern, the occupation of less- 
favored bridge sites is compensated by reduced adsorbate-adsorbate through-space repulsion 
and absence of metal sharing. Conversely, occupying the preferred hollow sites in the hf 
arrangements comes at the expense of metal sharing, which we have discussed as being 
unfavorable for Pd. 

On Rh (111) we find bb3, tb, and hf overlayer patterns to have the same -E c hem within 
0.03 eV. At low coverage, we find that Rh has a relatively weak inherent site preference, 
preferring t over b by 0.09 eV, and t over h and f by only 0.03 eV. Since both site preference 
and repulsion due to metal sharing are weak on Rh, we can explain the stronger -E c h em of 
bb3, tb, and hf as being due to not having any nearest-neighbor CO distances. The other 
four overlayers studied all have some nearest neighbors. 

On the (100) surfaces of Pt, Pd, and Rh, mixed site occupation is not observed at 
O=0.5. The square lattice geometry of the (100) surface permits single-site occupancy 
at O=0.5 without nearest neighbors. As a result (contrary to the (111) surface), a tb 
arrangement on the (100) surface would introduce adsorbates at shorter separations than 
in the tt2 arrangement. Therefore, on (100) surfaces we would expect overlayer patterns 
that minimize through-space repulsion and free of bonding competition most likely to be 
observed. Both tt2 and bb3 minimize adsorbate-adsorbate through-space repulsions and do 
not involve bonding competition, and this explains why we find E- m t to be no greater than 
0.06 eV for these patterns on all three metals. 
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V. CONCLUSION 



We have presented for the first time results for CO chemisorption at different coverages 
using DFT-GGA plus our first-principles extrapolation procedure for accurate chemisorp- 
tion energies. Our values for -E^em at higher coverage with CO occupying different sites and 
in different patterns show the experimentally observed structure to have the most favor- 
able -Echem within the selection of adsorbate arrangements modeled. We identify and dicuss 
adsorbate through-space repulsion, bonding competition, and substrate-mediated electron 
derealization as key factors determining preferred chemisorption patterns for different metal 
surfaces and adsorbate coverages. We rationalize how these effects, along with the inherent 
site preference energy for CO on each metal, balance to cause different trends in chemisorp- 
tion behavior as a function of CO coverage O on different transition metal surfaces. 
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